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ABSTRACT 

An  arbitrarily  orientable  Yagi-Uda  array  antenna  was 
modeled,  and  a  computer  simulation  run  to  obtain  the  input 
impedance,  gain  pattern  and  front-to-back  ratio  of  various 
arrays.   The  model  made  provisions  for  the  antenna  to  be 
operated  over  either  a  lossy  ground  plane  or  aboard  a  ship 
in  seas  of  specified  state.   Quick  solution  turn-around, 
with  CRT  display,  enabled  relatively  rapid  optimization  of 
numerous  arrays. 

Theory,  resultant  optimal  designs  and  performances, 
photographs,  and  program  listing  are  included. 
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I.   INTRODUCTION 

A.   BACKGROUND  OF  THE  STUDY 

The  most  concise  equation  that  describes  the  radiated 
electric  field  of  an  antenna  assumes  that  the  antenna 
operates  either  in  free  space  or  over  a  perfectly  conducting 
ground  plane.   Equations  in  this  category  are  simple  ones 
that  have  been  typically  used  in  text  books,  and  voluminous 
experimental  data  have  been  obtained  through  their  use. 
Emphasis  has  been  upon  simplification  of  equations  so  that 
problems  could  be  solved  manually  with  a  minimum  of  rigor. 
Equations  such  as  these  are  admittedly  inaccurate  because 
they  omit  the  component  of  radiation  produced  by  ground 
reflection  or  if  not  omitted  the  ground  surface  is  considered 
to  be  a  perfect  reflector. 

In  reality  the  simplest  antenna  operates  over  a  lossy 
ground  and  established  at  least  a  direct  and  a  ground 
reflected  wave.   When  simple  free  space  equations  are  modi- 
fied for  two  path  propagation  over  a  lossy  ground  plane  and 


^■The  ground  wave  is  composed  of  a  space  component  and 
a  surface  component.   A  particularly  important  point  may  be 
made  that  for  the  horizontally  polarized  wave  in  the  VHF/UHF 
range  the  strength  of  the  space  wave  greatly  exceeds  that  of 
the  surface  wave  so  as  to  render  the  surface  wave  negligible. 
This  then  allows  an  accurate  field  strength  to  be  calculated 
using  equations  which  assume  only  a  two-path  propagation 
model — the  direct  wave  and  the  lossy  ground-reflected  wave. 
Henry  R.  Reed  and  Carl  M.  Russell,  UHF  Propagation,  (New 
York:   John  Wiley  &  Sons,  Inc.,  1953) ,  p.  174 . 
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solutions  are  obtained  by  computer  with  a  Cathode  Ray  Tube 
(CRT)  display  of  the  gain  patterns,  interesting  results 
accrue. 

B.  STATEMENT  OF  THE  PROBLEM 

The  initial  problem  is  to  obtain  and  verify  the  equations 
required  to  describe  the  electric  field  of  a  multi-element 
linear  array,  particularly  of  the  Yagi-Uda  type.   A  recent 
work  which  provided  graphic  computer  solutions  to  single 
element  antennas  is  the  basis  of  this  study,  and  as  such 
this  study  is  a  follow-on. 2   Therefore  effort  centers  around 
adapting  the  existing  program  to  meet  the  requirements  of 
both  antenna  types. 

C.  OBJECTIVES 

The  objective  of  this  study  is  to  provide  a  near  real- 
time computer  graphic  solution  of  the  gain  pattern  of  an 
arbitrarily  oriented  Yagi-Uda  array  which  is  centered  above 
two  types  of  planes  (with  specified  e  &  a) : 

1.  the  lossy  ground  plane,  and 

2.  the  lossy  ocean  plane  which  rolls  and  pitches  the 
antenna  as  specified  by  the  sea  state. 

Examples  of  gain  patterns  and  other  output  parameters  are 
illustrated. 


2R.  W.  Adler  and  C.  B.  Robbins,  "The  Solution  and  Graphic 
Display  of  Gain  and  Patterns  for  Wire  and  Linear  Antennas  in 
the  Presence  of  Lossy  Ground",  Electrical  Engineering  Depart- 
ment, Naval  Postgraduate  School,  to  be  published  . 
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D.   SCOPE  AND  LIMITATIONS 

The  following  assumptions  have  been  made  for  the  study 
of  the  Yagi-Uda  array: 

1.  Propagation  is  confined  to  two  paths. 

2.  The  ship-ocean  model  does  not  make  any  provision 
to  augment  the  number  of  wave  paths  that  turbulent 
seas  might  produce. 

3.  The  current  on  the  elements  is  distributed  sinusoi- 
dally,  e.g.  the  elements  are  thin  (d  <_  A/100)  . 

4.  The  elements  may  be  spaced  arbitrarily  with  arbitrary 
lengths.   Thickness  can  be  changed.   Assignment  of 
the  element  lengths  is  constrained  by  the  fact  that 
if  elements  are  one  wavelength  then  the  solutions 
become  indeterminant. 

For  purposes  of  testing  the  resulting  equations,  various 
array  designs  were  attempted.   The  dimensions  of  the  array 
were  varied  using  the  method  of  iterative  search  to  uncover 
optimal  horizontal  and  vertical  designs  of  arrays  with  and 
without  a  reflector. 

Optimality  is  determined  by  three  criteria: 

1.  Input  Impedance  (Z^n) 

Where  it  was  possible  to  do  so,  the  array  was 
designed  to  have  a  reasonably  high  resistance,  e.g. 
£.  20ft  and  a  reasonably  low  reactance,  e.g.  _f_  10ft. 

2.  Front  to  Back  Ratio  (FBR) : 

Within  a  satisfactory  impedance  range,  the  FBR  is 
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.  .     3     . 

maximized.    This  expression  differs  in  some  cases 

from  the  expression  for  FBR  that  is  typically  used. 
3.   Power  Gain  (G) : 

Finally,  within  the  maximum  FBR  the  gain  is 
maximized. * 
Photographs  are  shown  of  the  linear  and  logarithmic 
results  that  are  obtained  from  a  variety  of  Yagi-Uda 
designs  placed  at  different  heights  above  ground  and  ocean 
environments . 

Preliminary  tests  were  made  which  showed  that  gain 
varied  with  change  in  element  thickness,  however  for  the 
results  obtained  throughout  this  study  the  element  thick- 
ness to  element  length  ratio  remained  fixed  at  1/200. 

Detailed  operating  procedures  are  found  in  Appendix  A. 


^For  purposes  of  the  study  the  following  expression 
for  FBR  isused 

FBR   =   G<9  ><t>)max    ,  /    °'l  4>    jl8  0* 

G(«^)max    /    1801  0    1360*       ■ 

This  differs  completely  in  some  cases  from  the  expression 
that  is  typically  used: 

FRR  -  G(e'90>   . 
FBR  "  G(e,270) 

4The  power  gain  used  is  measured  with  respect  to  an 
isotropic  source. 
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II.   THEORY  OF  THE  MODEL 

A.   ARRAY  IMPEDANCE  AND  ELEMENT  CURRENT 

The  radiation  problem  requires  the  solution  of  the 
individual  element  currents,  which  will  become  the  basis 
for  calculating  the  gain  of  the  array.   The  assumptions 
are  that  the  driven  element  is  excited  by  one  volt,  and 
that  the  antenna  is  thin  and  therefore  the  element  currents 
are  sinusoidal.   The  equations  for  the  current  column 
matrix  and  for  Imax  j.   are 


[I]   =   [Z]"   [V]  (1) 


Imax  i  =     .  Zi  (2) 

sin  b  a i 

[z]   is  the  combination  of  the  free  space  impedance  and 
image  impedance  which  is  found  by 

[z]   =     [Z]    +    C  [Z«]  (3) 

[z]   and  [z1]   are  NxN  dimensional — [Z']  being  the  mutual 
impedance  matrix  relating  the  actual  and  image  elements  as 
indicated  in  figure  1.   It  should  be  noted  that  the  matrices 
are  complex.   Complex  inversion  presents  no  particular 
problem  except  that  the  [z]  matrix  must  be  arrayed  according 
to  the  basic  arithmetic  operations  that  follow: 
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Actual 
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Image 
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>  y 


Figure  I.  Side  view  of  a  3  element 
Yagi-Uda  array  tilted  o* 
degrees  abovethe  ground. 
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<[R]  +  jCx])    •    ([ip]   +  j[iT])     =     [v] 


■R- 


(4) 


which  can  be  written  as  two  equations 


[R][iR]  -  CxlCij]  =  [V] 

[X][IR]  +  [R][IX]   =   [V] 


Re 


Im 


(5) 


which  can  be  solved  according  to  equation  (1) ,  where 
equation  (3)  now  appears  as: 


[Z]   = 


[R]   [-x] 

[X]    [R] 


(6) 


[z]  is  a  2Nx2N  matrix.   An  array  of  maximum  size  (5 
elements)  requires  inversion  time  on  the  SDS-9300  computer 
of  5  seconds—roughly  one-eighth  the  total  solution  time 
needed  for  the  complete  aximuth  and  elevation  patterns. 

The  reflection  factor  C  in  equation  (3)  is  a  function 
of  the  angle  which  the  antenna  element  makes  with  the 
ground.   C  is  written  as 


-iA» 


C   =   e~3   (Rh'cos  A'  +  jRv'sin  A') 


(7) 


The  reflection  coefficients  Rh '  and  R^'  are  the  values 
obtained  when  6=0  is  substituted  into  equation  (8) .   Theta 


M.  T.  Ma  and  L.  C.  Walters,  Power  Gains  for  Antennas 
over  Lossy  Plane  Ground,  ESSA  Technical  Report.  ERL  104-ITS 
74  (U.  S.  Government  Printing  Office,  Washington,  D.C., 
1969) . 
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is  equal  to  zero  because  the  coupling  which  takes  place 
between  the  free-space  and  the  image  element  (thereby 
producing  a  mutual  impedance  Z  * )  occurs  with  the  image 
directly  beneath  the  actual  element.   Equations  (3)  and  (7) 
show  that  the  horizontally  polarized  field  (A'=0  )  and  the 
vertically  polarized  field  (A '=90  )  give  weight  to  the 
value  of  [z1]  by  the  values  C=R^'  and  C=R  '  respectively. 
Since  the  solution  to  equation  (3)  is  actually  independent 
of  0  and  <J> ,  one  solution  for  this  equation  satisfies  the 
gain  expression  at  any  position  of  observation. 

The  values  of  R^  and  R^  which  are  also  used  in  the 
gain  equation,  are  a  function  of  the  observation  angle  6 
as  follows: 


Rh   =   (8a) 


cos  e 

- 

"2 

6 

A 

cos  e 

+ 

62 
6 

A 

cos   6 

- 

8 

A 

where 


Ry   =   ^ (8b) 

COS  9    +   75—  A 

*2 


2  * 

A  =   [1  -  (  L  sin  e)  ]  (9a) 

62 
*2   =   6[er  -  j  ±-f  Ob) 

W£q 

or  4  ** 

■     62   =   B[er  -  j  1'8    °    10  ]  (9c) 

r        f 

MH 


6 


Ibid. 
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The  matrix  [z ' ]  is  identical  to  equation  (6)  except 
that  it  is  written  in  terms  of  [R']  and  [X'].  The  geometric 
orientation  of  the  antenna  elements  are  central  to  the 
solution  for  the  mutual  impedance  matrix  [Z'].   Figure  2 
shows  an  arbitrarily  oriented  single  element  and  its 
corresponding  image  in  the  ground  plane.   Although  the 
true  coordinates  of  the  antenna  are  in  the  xyz  coordinate 
system,  the  equation  for  solving  for  the  impedance  fixes 
the  primary  element  along  a  vertical  axis — here  it  is  shown 
to  be  z11.   The  secondary  element  becomes  the  image  element 
and  is  separated  from  the  primary  element  by  the  distances 
ZQ  and  Y  .   The  linkage  between  the  two  coordinate  systems 
xyz  and  z'y'z1  is  01.   The  equations  for  ZQ  and  Y   are 
observed  to  be 

(10a) 
(10b) 
(11a) 
S  ,  =   -S  sin  29'   =   -S  sin  2A'  (lib) 

y 

The  solution  of  self  impedance  is  performed  by  assuming 
that  Z  =0  and  YQ=  2  r,  where  the  radius  r  is  assumed  equal 
to  &j/200.   (This  makes  Jt-/d  =  100,  which  for  i±   =  A/4 
satisfies  the  thin  antenna  specifications  of  d  <_  A/100) . 
The  mutual  impedances  of  the  free-space  array,  where  the 
elements  are  all  parallel,  are  solved  using  Z  =0  and  YQ=y^ , 
the  separation  between  the  two  elements  in  question.   Under 
dynamic  ship  motion,  which  affects  6'  and  h,  the  antenna 
height  above  the  earth  plane,  the  mutual  impedance 


zo 

= 

-2h  cos  9 ' 

= 

-2h  sin  A1 

Yo 

= 

2h  sin  9 ' 

= 

2h  cos  A1 

sz' 

= 

S  cos  26 ' 

= 

S  cos  2A' 
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actual 
(primary) 

antenna   L 


image 
(secondary) 
antenna 


Figure  2.  The  geometry  of  an  arbitrarily  oriented  element 
in  free  space  and  its  image  (used  for  solution 
of  mutual  impedancesZ1) 
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matrix  [z1]  will  be  affected  according  to  equations  (10)  and 
(11) ,  but  the  self  and  mutual  impedances  of  the  [z]  matrix 
will  not  be  affected  since  it  represents  the  free-space 
array  impedance. 

Impedances  are  solved  according  to  the  following: 
Ls/2 


R  =  -30 


ZlFf 


^l(z0+sz.+ip)+^^^(^sz'-^ 


-V2 

-2  8in(2irr)  ^^       J 
2X  J 


7\ 


z  _0 

2X 


Y+Sy' 


2  sin(2rrr)  gog/^to)  -  sin  (2^)  _  sin  (2^) 


•  V 


sin  (2v JL  -  I si) 
2\ 


ds 


(12a) 


Ls/2 
X  =  -30    f 

-Ls/2 


cos(2uri)   (Z„+Sz'+^) 
rl 


cos(2iTr2)        ,  L^x 

+  =-  (Z0+SZ'-_P)  - 

r2       °      2X 


2  cos(27Tr) 


cost £) 


(z0+sz') 


(Y0+Sy') 


-cos(2TTr^)    cos(2iTr2) 


2  cos(2TTr)     ,ttL~.. 
cos  ( E) 


LS 


sin  (2iT-2y  ~  I  s  I  ) 


ds 
(12b) 


'H.  C.  Baker,  and  A.  H.  Lagrone;  "Digital  Computation 
of  the  Mutuai  Impedance  Between  Thin  Dipoles;"  IRE 
Transactions  on  Antennas  and  Propagation;  March,  1962; 
AP-10,  No.  2;  pps  172-178. 
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P2   =   ^Yo+Sy')2  <13a> 

r   =   [p2  +  (Z0+Sz«)  ]  (13b) 

r  2  Ln  2% 

rl   =   EP   +  (Zo+Sz'+2^)  ]  (13c) 

r0   =   [p2  +  (Z  +S  »-5R)2]  (13d) 
2              o   z   2 


Where  the  distances  r,  r..  ,  and  r  ,  are  the  respective 
distances  in  figure  2  from  the  bottom,  center  and  top  of 
the  primary  antenna  to  the  differential,  ds,  on  the 
secondary  antenna. 

B.   THE  FAR  FIELD 

To  solve  for  the  electric  field,  a  point  of  observation 
is  chosen,  and  each  element  in  the  array  is  viewed.   Each 
element  contributes  to  the  total  field  according  to  its 
vector  phase  and  amplitude,  its  separation  from  some 
reference  element,  the  tilt  angle  of  the  array,  and  the 
bearing  to  the  point  of  observation.   Figure  3  shows  the 
geometry  of  a  2  element  horizontal  Yagi-Uda  tilted  a  degrees. 

The  equation  of  the  field  produced  by  one  element  located 
at  the  origin  is  first  developed.   Then,  the  presence  of 
other  elements  in  the  array  is  taken  into  account. 

The  equation  at  the  point  P  for  an  element  at  the  origin 

jU  (t    Q> 

Ti     =   V-  sin[6(£i±x)]  e  (14) 

Current  flowing  in  the  incremental  length  dx  produces  a 
field,  and  there  will  be  a  phase  associated  with  it 
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ith 
element 


Array  axis  (ZY  Plane) 


reflector 
element 


P'atoo 


Figure  3.  Yagi-Uda  array  with  horizontal  elements  and 
a  tilted  axis,  as  viewed  from  a  distant  point  P. 

(ay=O,0'=9O) 
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with   respect   to   point   P   as    the   differential    length  moves 
from      -—     to      +   —5.    .      The   delay    is    x    cos    \\> ,    which    is    found 
by   the   dot   product   of   OX   and   OR,    as    follows: 


OA    =     OX  •  OR 
where. 

^  A  *  •* 

OX    =     x  D-  si*1  9 '   cos  *  '  +  J  s^cl  6  '   si*1  <J> '  +  k  cos  6  ']  (  15a) 

^  *  *»  ». 

OR    =        [i  sin  e     cos  <f>    +  j  sin  6     sin  <j>     +  k  cos  9  ]  (  15b) 

so   that 

OA  =  x  [sin  6  sin  e'  cos  (<t>-<f>')  +  cos  9  cos  9']     (16) 

or 

OA   =   x  cos  if)  (17) 

where 

cos  ip   =   sin  9  sin  9'  cos  (^-cf)')  +  cosS  cos  91  •  (18) 

It  is  now  clear  that   s  =  r  -  x  cos  i>      and 

■inT*.  t^~X    cos  i>    \  -i 
I.   =   IQ   sin[6(£.±x)]  eDw[t-(— H }  ]  <") 

where 


(  i  +   X   ,    X  <  o 
±x   =  { 

(  t  -  X   ,    X  >  o 


Although  there  is  a  phase  difference  between  r  and  s,  at 
points  taken  far  from  the  origin  their  magnitudes  are 
approximately  equal. 

Observe  that  by  translating  the  xyz  origin  as  well  as 
the  off-axis  element  to  the  point  B  located  along  the 
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observation  vector  r,  the  magnitude  of  the  new  observation 
vector  relative  to  the  new  origin  is  BR,  which  differs  in 
length  from  the  old  observation  vector  OR.   The  amount  of 
change  equals  that  which  affects  the  observation  vector  s. 
Furthermore  the  angle  by  definition  remains  unchanged. 
Therefore,  the  expression   r  -  s  =  x  cos  i>      holds  for  every 
radiating  element  in  the  array,  not  just  for  the  element  at 
the  origin.   From  this  it  can  be  concluded  that  the  integral 
equation  for  the  field  is  identical  for  each  element  in  the 
array. 

q 

Following  the  method  used  by  Kraus: 

I. 

„                  j    60    tt    sin    <J>  f2-  ,  ,„,  , 

E<t,i      =      J /  I±   dx         .  (21) 


sX 


1 


O  TT 

Substituting  for  I  and  letting   6  =* —  =  _   ,  we  see  that: 

X     c 

I. 

j  30  8  Ju(t-rr)    7  •  zoo  ±o  s    jSxcos^  ,   ,00v 

E<j».  =  ^ sin  4»  I_  e     c     /  sin(B£.±8x)eJ      dx.(22) 

1  r  i  -l  J 

When   this   expression   is    integrated   the   result   is: 

j  60  sin  <f>  I0.   eju(t"^ 
E*i     ~    h [cos(6£-cos  ij;) -cos  ( 6^)1      .    (23) 

•    2    , 

r  sin     i> 


8J.D.  Kraus;  Antennas;  (New  York:  McGraw-Hill  Co.  1950), 
pp.  135-141. 
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For  an  array  of  elements  the  field  is  the  superposition  of 
the  fields  of  each  element  taken  separately.   E.  is  then 
the  sum  of  E<{>.  given  by 

N 


e"3    r    60    sin    4>V*I^      Ccos  (B£  •  cos^)  -cos  (B4 _• )  ]  .     (24) 

E<j>     "  ~T~,        Z-j     i 

*  r   sin      <{»        T~ 

i=l 

It  is  recognized  that  the  time  term,  e~^    has  been  dropped 
and  that  B  has  been  substituted  for  u/c. 

The  array  factor,  or  the  phase  delay  OB  introduced  by 
the  separation  between  elements,  is  next  taken  into  account. 
From  figure  3,  using  the  reflector  as  reference,  the 
separation  between  elements  as  viewed  from  point  P  is 
OB  =  OY'OR  ,  where 

OY   =   jy-cos  a  +  ky-sin  a  (25) 

A  A  A.  * 

and  OR  =   i  sin  6  cos<t>  +  j  sin  esin<t>  +  k  cose  (26) 

gives         OB   =   y^  (cosa  sin  6  sin  <J>  +  sin  a  cos  6  )  .  (27) 

If   cos  fi  =   cos  a  sin  6  sin  <J>  +  sin  a  cos  8  (28) 

then  OB   =   y.  cos  n  (29) 

represents  the  phase  lead  of  the  i"*-*1  element. 

Having  found  the  slant  separation  between  elements  to 
be  expressed  by  equation  (29) ,  the  current  for  any  element 
in  the  array  may  be  written  as 

I     =   i      eJBYiCOs  n  a  (30) 

o^      Max  i 

The  phasor  sum  of  currents  in  the  array  becomes 
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N 

-  E 

i=l 


jBy^COS  Q 
I0   =  2-/    Imax  i  e  •  (3D 


Upon  substitution  of  equation  (31)  into  equation  (24)  the 
result  is 

-ifir  N 

j  60  e  J   sin  j  V^  T       JByicos  n 

.  E<j>  .      ~                         2  /  ->      max    i 

1                    r   sin      t|>  i=l 


[cos(6ai.Cos    <M    -   cos(8£i)] (32) 

where  Imax  a    is  found  from  equation  (2) . 

Equation  (32)  represents  the  free  space  radiation  of 
the  horizontally  polarized  Yagi-Uda.   The  equation  does  not 
consider  the  ground  reflections.   Figure  4  shows  the  direct 
and  ground  reflected  waves,  and  the  geometry  shows  the  dif- 
ference between  their  path  lengths.   The  element  separation, 
written  as: 

2h±   =   2(h  +  y.  sin  a)  (33) 

Ai 

is  a  function  of  the  path  difference  vector  AB,  since 
AB  =  2hicos  6  . 


Thus 


E^         =      E^    +   Rh  E+    e-326hic°s   9   .  (34) 


And  finally,  equations  (32)  and  (34)  are  combined  to  produce 
the  final  space  wave  in  the  phi  direction: 

f    =   j  60  e  J    sin  4) 

T        ■       2 

r  sin  \\j 


:  [1  +  Rh   e-j2Bhicose]- 

max.;    L  h  J 

jByicos  ^cos(6£iCOS^)    -cos(3£i)].     (35) 


•i 
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actual  antenna 
(end  on) 


image  antenna 
(end  on) 


Figure  4.  The  tilted  Yagi-Uda  receiving  horizontally 

polarized  waves  via  direct  and  reflected  paths 

(The  array  is  tilted  to  present  less  confusion 
in  distinguishing  angles) 
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Observe  that  when   cos  fi  =  sin  0  sin  4>  and  h^=h  equation 

(35)  agrees  with  Ma  and  Walters  equation  for  the  horizontal 

9 

Yagi-Uda.    The  equation  for  the  0  component  of  the  horizon- 
tal Yagi-Uda,  Ee,  is  not  derived  with  the  same  rigor  since 
it  represents  the  cross  polarized  radiation,  which  exists 
only  off -axis  and  is  of  minor  interest.   The  equation  is: 

F    _   j  60  e  J   cos  0  cos  j> 

^T  .  2  , 

r  sin   ^ 

^A        r        -j26hiCOS  0_ 

>    I      [1  -  R  e  J    1      ] 

/  j       max^       v 

eJ6yiCos  fi[cos(6£.cos  *)  -  cosCB^)]  .      (36) 

10 

The  equation  for  the  vertical  dipole  is  well  known. 

If  R   is  used  insteady  of  R,  then  equations  (31)  and  (34) 
may  be  applied  to  the  standard  free  space  equation  to  arrive 
at  the  general  solution  for  the  vertical  array: 


0T  JL-j     max^ 

r  sin  o    f=7 


JBYicos  ^[cos(fU.cos  „,)  _  cos(0£i)]        (37) 


Specifically/   cos  i>    =  sin  0  sin  4>   and   cos  ^  =  cos 
since   a  =  0  =  0   for  the  vertical  Yagi-Uda. 


g 
M.T.  Ma  and  L.C.  Walters,  op.cit .  p.  41. 

l°Kraus,  loc.cit . 
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C.   THE  SHIP-OCEAN  MODEL 

Figure  5  shows  the  effects  of  the  dynamic  ship  model 
upon  the  Yagi-Uda.   Changes  in  the  sinusoidal  waves, 
considering  that  a  roll  of  eight  degrees  per  sea-state  is 
produced  and  a  bow  pitch  of  three  degrees  per  sea-state  is 
produced  (representative  of  a  light  cruiser) ,  will  affect 
the  angles  e'  and  a  as  follows: 

e1   =   et'  -  (a-j^cos  4)'  +  A2sin  <j>  ' )         (38) 
a        =      at   -    A]_sin  <j>  '  +  a2cos  <|>  *  (39) 

where 

A,   =   (wave)  sin  (course)  (40a) 

A2   =   (wave). 3  cos  (course)  (40b) 

and  where 

wave   =   8 (sea  state)  sin  (y)   .  (40c) 

The  effective  height  of  the  array  will  change  according 
to 

h  =   h.  cos  A-l  cos  A-  (41) 

It  should  be  noticed  that  the  parameter  <f>  '  until  now 
has  been  somewhat  extraneous  for  the  solution  to  antennas 
over  land,  since  variation  of  the  observation  angle  <£ '  and 
the  angle  <j>  accomplish  the  same  thing.   However  the  angle 
<J>  '  now  is  significant  in  specifying  the  antenna  orientation 
with  respect  to  the  ship  heading,  which  then  determines  the 
aspect  that  the  wave  presents  to  the  antenna. 
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New  j 

Array        r 


Ships 
y     Heading 


Figure  5.  Orientation  of  the  Yagi-Uda 

aboard  a  ship.  Angles  a,  and  az 
represent  roll  and  bow  pitch 
respectively. 
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D .   GAIN 

The  gain  of  an  antenna  is  equal  to  the  ratio  of  the 
power  intensity  to  the  power  density,  and  is  expressed  as 
follows: 

G  =  i^r  (42) 

Win 


where 


or 


Wr   =   rz  (  E  x  H*)   watts/solid  angle 

2       2 

W    =   JE |E|  (43a) 

r      120* 


2 
in      o'    in 


and  W4_  =   |Ib|   R._   watts  .  (43b) 


The  final  expression  for  gain  appears  as 


G  =  r2lEeTr  +  Ie^tT  (44) 

2 


30  |Ib|   R 


in 


The  equation  for  gain  becomes  independent  of  r  when  the 
equations  (35)  and  (36)  are  substituted  into  the  equation 
above . 
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III.   RESULTS  OF  ARRAY  PERFORMANCE 

A.  RESULTS  OVERVIEW 

Basically,  four  types  of  arrays  were  examined  for 
optimization:   horizontal  and  vertical  with  and  without 
reflectors.   Experiments  were  conducted  to  determine  the 
design  of  arrays  that  afford  optimal  performance  as  defined 
in  the  Scope  and  Limitations  section  of  the  study.   Starting 
with  the  two  element  array,  the  parasitic  element  function- 
ing both  as  a  director  and  as  a  reflector,  the  element 
spacings  and  lengths  were  varied  to  obtain  optimization. 
The  arrays  were  lengthened  by  adding  one  element  at  a  time, 
up  to  a  total  of  five  elements,  each  time  solving  for  an 
optimal  design.   Numerous  iterations  were  required,  and  the 
time  required  to  optimize  an  array  varied  from  ten  minutes 
for  a  two  element  array  to  slightly  over  three  hours  for  a 
five  element  array.   (The  vertical  array  was  more  difficult 
than  the  horizontal,  and  the  array  without  reflector  was 
more  difficult  than  that  with  reflector.)   Not  all  of  the 
multi-element  arrays  were  optimized  according  to  the  rigid 
criteria  chosen  for  this  study.   Specifically,  the  vertical 
designs  having  four  and  five  elements  were  arrived  at  by 
relaxation  of  the  impedance  threshold  criterion  thereby 
producing  a  suboptimal  design. 

B.  THE  TWO  ELEMENT  ARRAY 

Referring  to  Table  I  for  the  horizontal  two-element 
array  with  a  director,  maximum  FBR  is  maintained  when,  for 
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TABLE  I 

OPTIMAL  2 -ELEMENT  HORIZONTAL  ARRAY 
AT  VARIOUS  ELEMENT  SPACINGS  USING 
THE  PARASITE  AS  A  REFLECTOR  AND  AS 
A  DIRECTOR 

(Gain  and  Front-to-Back  ratio  in  db) 
(Spacing  and  length  in  M 


d2               .012              .100         .150         .200         .251         .295 

D 
I 
R 
E 
C 
T 
0 
R 

GP 

FBRp 

di 

i                       i                   J                  i                  i 

10.3        .      10.6       JlO.l        i    9.5           8.7       !    7.8 

13.3       i      11.9       »    6.7        ;    3.7       j    1.9       J    0.8 
.492,            .460  J       .449'       .438!       .424]       .390 

R 
E 
F 
L 
E 
C 
T 
0 
R 

G 
P 

FBRp 

10.5       j      10.9       '10.7        110.5       -10.1           9.9 
10.9       •      13.6       j 13.6       |l3.3       J12.5       Jll.9 
. 501  |            .50l|       .50l!       .50l!       .506|        .501 

f   =146.   h=lm.   In  both  cases  the  length  of  the  driven 

MH  * 

element  is  A/2.   (L^  and  L   are  director  and  reflector 

lengths  respectively.)   (The  subscript  p  denotes  peak, 
e.g.  the  length  of  the  parasitic  element  was  varied  until 
the  peak  FBR  was  first  obtained,  then  within  that  value, 
the  gain  peak  was  found  by  making  finer  adjustments  to 
length. ) 
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increasing  element  spacing,  there  is  a  corresponding  decrease 
in  director  length.   Apparently/  best  performance  occurs  when 
spacing  is  .  01X  (.021m)  because  this  is  where  the  maximum 
FBR  appears.   However,  no  mention  has  been  made  about  the 
affect  upon  impedance  of  varying  the  element  spacing. 

In  general  it  was  found  that  the  resistance  and  inductive 
reactance  decreased  together  (from  70  to  40  ohms  and  from 
40  to  10  ohms  respectively)  as  elements  were  placed  closer 
together.   When  the  spacing  closed  to  .  05X,  the  currents  in 
the  elements  evidently  became  very  large,  and  the  value  of 
resistance  became  less  than  one  ohm  (too  low  for  convenient 
matching) .   Thus  the  FBR  and  G  are  insufficient  criteria  of 
optimality  since  impedance  can  fall  (in  some  cases  may  rise) 
to  unreasonable  values.   Optimality  in  the  case  of  the  di- 
rector array  in  Table  I  actually  occurs  at  .IX  (.205m) 
spacing . 

When  the  parasitic  element  is  used  as  a  reflector  it  is 
not  necessary  to  change  the  reflector  element  length  as  the 
spacing  between  elements  is  changed  in  search  for  the 
maximum  FBR.   The  maximum  FBR  occurs  at  . 5X  (1.03m)  for  all 
spacings  except  at  .25X  (,515m)  j  that  is  , the  maximum  FBR 
occurred  when  both  element  lengths  were  the  same. 

The  gain  curves  shown  in  figure  6  differ  from  other 
published  experimental  and  theoretical  data  (and  there  are 
fairly  wide  variations  in  these  results  too) ,  inasmuch  as 
the  value  of  the  peak  gain  is  somewhat  higher  than  that 
typically  recorded.   This  is  because  published  data  typically 
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14 


Element  Spacing 
i i 


.1  .2  .3 

Figure  6.  Comparison  of  gain  vs 
element  spacing  for  optimal  2— 
element  horizontal  arrays. These 
curves  represent  data  in  table  I, 
where  the  array  is  a  reflector- 
driven  element  combination  in  one 
case  and  a  director-drivenelement 
combination  inthe  other  case. 
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does  not  reflect  a  lossy  ground  plane  (two-path  propagation 
is  ignored).   Also,  the  gain  of  the  reflector  array  exceeds 
that  of  the  director  array,  which  generally  disagrees  with 
other  findings.   Again,  the  explanation  for  these  differences 
is  tracable  to  the  assumption  that  a  lossy  ground  wave  is 
being  propagated  along  with  the  direct  wave. 

C.   THE  MULTI-ELEMENT  VERTICAL  AND  HORIZONTAL  ARRAYS 

Tables  II  and  III  list  the  optimal  and  suboptimal  design 
parameters  for  the  fm,=3  0  arrays.   Notice  that  the  four  and 
five  element  vertical  arrays  have  relatively  wide  reflector 
spacings--these  are  suboptimal  designs.   By  comparison  with 
the  optimal  design  that  was  found  for  the  four  element 
vertical  array  (where  R^n  >_  20ft)  the  FBR  is  16.4  db  whereas 

the  suboptimal  design  (Rin  <  20ft)  produced  a  FBR  of  22.4  db 

12 
—  quite  a  drop.   The  two  designs  differed  completely. 

Furthermore  they  responded  to  slight  parameter  variations 

differently. 

The  four  element  vertical  array  was  fine-tuned  for 

maximum  FBR  only  to  discover  that  the  impedance  was 

unacceptable;  all  subsequent  fine  adjustments,  introduced 

for  purposes  of  raising  the  resistance  to  a  value  above  20ft, 


11C.  R.  Fry,  The  Yagi-Uda  Aerial — A  Short  Design  Review 
and  Bibliography,  (Valcartier ,  Quebec :  Canadian  Armament 
Research  and  Development  Establishment,  May  1966),  p.  23. 

l2The  optimal  design,  as  compared  with  the  suboptimal 
design  shown  in  Table  II,  was  Lj_=  .49  ,  .  47  ,  .  45  ,  .  43  ;  dj_=.22, 
.19, .15  (L  and  d  in  X) . 
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caused  a  drastic  reduction  in  FBR — below  16  db.   Thus,  this 
particular  design  was  abandoned;  another  completely  different 
design  provided  an  acceptable  impedance  with  a  peak  fine- 
tuned  FBR  value  of  16.4  db  (compare  parameters  in  footnote 
#12  with  the  comparable  parameters  in  Table  II) . 

Except  for  the  two  larger  vertical  arrays  (four  and  five 
elements) ,  all  changes  that  were  made  for  the  purpose  of 
increasing  R^n  or  reducing  X^n  inevitably  improved  the  FBR 
figure  and  in  none  of  these  cases  did  the  manipulation  of 
parameters  for  purposes  of  adhering  to  the  impedance 
criterion  result  in  obtaining  a  lower  FBR  or  G  value. 
Because  such  manipulation  did,  however,  cause  a  drop  in  the 
FBR  for  the  two  vertical  arrays  they  are  considered 

anomalous . 

13 
The  optimal  horizontal  and  vertical  arrays   have 

tapered  elements  in  the  region  next  to  the  driven  element, 

but  the  last  few  end-directors  of  the  four  and  five  element 

horizontal  array  show  a  slight  inverse  taper.   Writers, 

however,  generally  agree  that  elements  and  spacings  that  are 

gradually  tapered  towards  the  end  of  an  array  will  usually 

give  best  results. 

It  is  obvious  that  the  optimum  parameter  measurements 

obtained  for  the  reflector-driver  combination  (two-element 

array)  do  not  necessarily  ensure  optimization  of  successive 

experiments  involving  the  addition  of  another  director.   A 


■^The  optimal  arrays  include  all  horizontal  plus  the  two 
and  three  element  vertical  arrays. 
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design  readjustment  of  all  parameters  is  required  each  time 
the  array  length  is  changed  (e.g.  each  time  an  element  is 
added  or  removed) .   A  corollary  to  this  for  arrays  with  more 
than  two  elements  is  that  the  FBR  will  decrease  with  the 
addition  of  a  director  whose  length  and  spacing  are  identical 

to  those  of  the  preceding  director.   This  was  observed  in 

14 
testing  the  uniform  array   at  fMH=10  where  the  FER  dropped 

from  14.6  db  as  a  four  element  array  to  10  db  as  a  five 

element  array. 

With  regard  to  director  spacings,  perhaps  it  is  more 
than  just  coincidental  that  the  designs  which  produced 
optimization  of  the  longer  horizontal  arrays  resulted  in 
equispaced  director  elements.   Nothing  could  be  obtained 
to  substantiate  whether  other  analysts  would  agree. 

From  Table  III,  comparison  between  the  horizontal  three 
and  four  element  arrays  as  well  as  comparison  between  the 
four  and  five  element  arrays  shows  that  the  FBR  improvement 
is  greater  than  the  G  improvement  as  array  lengths  increase. 
An  added  advantage,  besides  increasing  FBR  and  G,  is  that 
the  greater  the  number  of  elements  used  in  the  optimal  array 
the  greater  the  tendency  for  automatic  reactance  cancella- 
tion.  While  the  reactive  components  of  the  two  and  three 
element  arrays  are  small  they  did  not  become  so  as  a 
coincidence  of  optimization;  manipulation  of  parameters  was 
necessary  to  obtain  reactance  decreases  in  conjunction  with 


14Design:  LL=.  533  ,  .  500  ,  .  434  ,  .  434  ,  .  434  and  dj_=  .  233  ,  .  100  , 
100,  @  h=3  (L,  d  and  h  are  in  X) . 
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the  need  to  increase  resistance  (and,  as  inferred  earlier, 
these  manipulations  in  the  case  of  the  two  and  three 
element  arrays  produced  slight  decreases  in  FBR  and  G   ) . 

Comparison  between  the  optimal  horizontal  array  and  its 
vertical  counterpart  indicates  that  the  horizontal  array 
consistently  out-performs  the  vertical  array  in  FBR,  but 
matches  it  in  G.   Furthermore,  when  the  optimally  designed 
horizontal  array  is  operated  in  the  vertical  position  and 
vice  versa  a  degradation  in  performance  results.   An  improve- 
ment of  almost  16  db  in  FBR  and  1.5  db  in  G  can  be  obtained 
when  the  array  is  optimized  in  the  position  for  which  it  is 
intended  to  be  operated. 

Table  IV  gives  the  results  of  operating  the  optimal 
array  without  its  reflector;  it  is  interesting  to  observe 
the  importance  that  the  reflector  has  in  determining  the 
overall  array  performance.   Practically  no  directivity  was 
obtained  without  a  reflector  even  with  a  four  element 
(three  director)  array.   The  unusually  large  G  that  results 
for  the  three  and  four  element  vertical  arrays  are  anomalies 


l^Fishenden  reiterates  the  popular  belief  that  a 
convenient  method  of  altering  the  input  impedance  without 
affecting  FBR  or  G  is  by  varying  the  spacing  between  the 
reflector  and  the  driven  element  from  A/8  to  A/4.   Such 
attempts  in  this  study  did  however  noticably  affect  the 
FBR  and  G  (particularly  in  the  case  of  the  vertical  four 
and  five  element  arrays,  which  were  particularly  trouble- 
some) .   R.  M.  Fishenden  and  E.  R.  Wiblin,  "Design  of  Yagi 
Aerials",  Proceedings  IEE,  96,  Pt.  3,  (Jan  1949) ,  p.  6. 


41 


TABLE  IV 

THIRTY  ARRAYS  DESIGNED  AS  SHOWN  IN  TABLE  II 
BUT  WITH  REFLECTORS  REMOVED 


Horiz  @  Horiz 
Orientation 

■ 

Vert  @  Vert 

Orientation 

N    1      2      3      4 

12     3     4 

G   [      7.8   11.5    5.0    7.9 
FBR  i      0      0      1.7    1.8 
Rin  I     60.6   29.4    4.1   33.0 
Xin  i    -12.3   -0.1-15.3   -9.0 

7.8    4.8   15.5   15.6 

0      0      2.8    2.9 

60.6   34.2   18.1   16.1 

-12.3   -9.7   -6.4   -4.2 
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and  cannot  be  explained  (undoubtedly  it  is  connected  with 
the  fact  that  these  two  arrays  are  actually  suboptimally 
designed  to  begin  with.) 

Such  poor  performance  is  not  without  remedy,  however, 
as  the  results  show  in  Tables  V  and  VI.   Arrays  without 
reflectors  may  also  be  optimized  with  considerable  improve- 
ment in  performance.   Unlike  the  arrays  with  reflectors, 
the  optimal  design  of  arrays  without  reflectors  gives 
identical  optimal  results  regardless  of  whether  the  array 
is  operated  in  the  horizontal  or  vertical  position.   Inas- 
much as  the  reactance  is  much  greater  than  10ft,  these 
arrays  do  not  represent  optimal  designs  according  to  the 
purposes  of  this  study,  yet  none  better  could  be  found. 

D.   SELECTED  INTERESTING  PATTERNS 

Numerous  photographs  were  taken  of  the  following  antenna 
patterns  (in  a  few  cases  comparisons  between  the  linear  and 
the  log  gain  plots  are  shown) : 

1.  a  three  element  array  operating  over  dry  land,  fresh 
water  and  sea  water, 

2.  a  three  element  array  at  different  heights  over  land, 

3.  several  three  element  arrays  operating  on  a  ship 
underway  in  heavy  seas, 

4.  a  three  element  array  with  its  axis  tilted  at  various 
angles  from  the  horizon  and 

5.  some  resulting  patterns  when  an  array  is  operated 
off-frequency . 
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TABLE  V 

PARAMETERS  OF  OPTIMAL  30  MHz  ARRAYS  WITHOUT  REFLECTOR 

(Horizontal  &  Vertical  Optimal  Arrays  Identical) 

(L  and  d  in  A) 


2-Element 

3-Element 

4-Element 

H& 

OV 

L 

.49    .48 

.51    .475    .475 

.52    .44    .45    .43 

RE 

i 

IR 

ZT 

d 

.05 

.08     .10 

.13    .13    .13 

TABLE  VI 

RESULTS  OF  OPTIMAL  ARRAYS  WITHOUT  REFLECTOR 

(Specifications  and  parameters  same  as  in  Table  V) 
(G  and  FBR  in  db,  R  and  X  in  ohms) 


N   2 


1 

G    »  19.7   13.6   10.1 
I 
i 

FBR  !   0     14.9   26.2 


R    •   4.8   76.6   73.9 
i 

i 

X    !  -6.2   57.4   74.5 
I 
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Figure  7   Optimal  30  MHz 
3-element  horizontal  array: 
linear  plot  (a)  and  log 
plot  (b) ,16 


Gp 
FBFL 


R 


in 


=  12.6  db 
=  22.0  db 
=  31.6  Jl 


Lm 


=  .l-n. 


Figure  8   Optimal  30  MHz 
4-element  horizontal  array: 
linear  plot  (a)  and  log 
plot  (b) .16 


F§Rj 
Rin 


13.3  db 
23.9  db 
36.1.n-  X 


in 


=  -2.3-0- 


Figure  9   Optimal  30  MHz 
5-element  horizontal  array: 
linear  plots  la)    &  (b) ,  and 
log  plot  (c) .I6 


Gp 

FBRp 

Rin 


-  13.5  db 
=  27.3  db 
=  28.3iV   X 


in 


=  -1.0-n. 


16 


Parameters  that  are  not  shown  in  Table  II  but  which 


are  common  to  all  the  30  MHz  arrays  are:   ht=3,  e=5.   a=10 
$  =90.   6  =85,  however  figure  9  (c)  differs  in  observation, 
with  <j>=20.  and  9=60. 


-3 
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(a) 
Dry  Land 
Gp    =  10.7  db 
FBR„  =   9.6  db 


(b) 
Fresh  Water 


F§RT 


=  9.9  db 
=  9.5  db 


(c) 
Sea  Water 
G     =  12.7  db 
FBRn  =   9.5  db 


Figure  10.   Ten  MH  three  element  vertical 
array  radiating  over  three  types  of  terrain. 


Figure  11.   150  MH  Optimal 
three  element  horizontal 
array — linear  plot  (a) ,  and 
log  plot  (b) .is 


GP 

™! 

Kin 

xin 


■■   12.9  db 
:  22.0  db 
=  28.7  Ji 
=-13.5  jt 


17Li=.506f  .500,  .450;  dj_=.  233  ,  .133  ;  ht=2  .  ;  fMH=10 . ; e=10 . , 
a=10"3f  4>p=90.,  ep=84.  (land)  ;  e=80.,  a  =  2 .  ,    <J>  =90.,  6p=76 . (H2O) ; 


e=80.,  a=10-2 


,  <f.p=90., 


j=75.(sea).  (Lf  d  ana  h^  are  in  X) 


18L.j  =  .510,  .490,  .460;  di=.250,.08;  ht=15;  e  =  5.;  a=10 
4>p=90.;  ep=89.  (L,  d  and  ht  are  in  X) 


=i n~3. 
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PS* 


=  80.0 


ht 
G 

FBRp 


(b) 
=  2.0 
=  10.7  <a 
=   9.6 


,=84.0 


,=76.0  G. 


FBRr 


(d) 
=  2.0 
=  13.3  @ 
=  27.3 


=  83 


Figure  12.   Ten  MHz  3-element  vertical  arrays  (a)  &  (b) 
and  horizontal  arrays  (c)  &  (d)  at  different  height 
over  land.19   (G  and  FBR  in  db) 


19 


See  figure  16  for  parameters.  (ht  is  in  X.) 
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(a) 

W 

(c) 

uo   =    0. 

u   =    30. 

(jo    =    90. 

stable 

pt.roll   =    15.4° 

st. roll   =    31.0* 

up  pitch   =      5.5* 

dn. pitch  =   11.0° 

e'   =  71.5 

6'    =    120.0 

a      =    64.0 

a      =    -11.0 

h      =      1.82 

h      =         1.68 

Figure  13.   Ten  MH  three  element  horizontal  array 
aboard  ship  in  rough  seas.   Sea  state  =  7  and  relative 
^■ir-or'-H  on  of  waves  =  40.   Gain  varies.  2^ 


(a) 

ID  J 

(C) 

co    =    0. 

oj    =    30. 

(jo    =    90. 

stable 

pt.roll   =    21.0° 

pt.roll   =   41.5° 

up   pitch  =      3.6° 

up. pitch  =      7.2* 

8 '    =    69.0 

6'    =    48.5 

a      =      3.6 

a      =      7.2 

h      =      2.80 

h      =      2.20 

Figure  14.   Thirty  MH  three  element  horizontal  array 
aboard  ship  in  rough  seas.   Sea  state  =  6  and  relative 
direction  of  waves  =  60.   Gain  varies. 21 


20Li=.506,  .500,  .450;  di=  .  233  ,  .  133  ;  e  =  80.;  o=5.;  4>  =  90.; 
=8  0.;  ht=2.  (L,  d  and  h-t  are  in  X) 

21L,=.51, .50, .46;  di=.25,.08;  e  =  80.;  a=5 .  ;  4>=90.; 
=85.;  nt=3.   (L,  d  and  h-t  are  in  X) 
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(a) 

(bj 

=  0. 

=  30. 

stable 

15.4*p.roll 
5.5  u. pitch 
9'  =  58.5 
a   =   5.5 
h   =  14.5 

Figure 

15. 

150 

MHz  3-eleme 

=  90. 

31. °p. roll  27."  s. roll 

11. °u. pitch  9.6°d.pitch 

e»  =  59.0  e-  =  i2o. 

a     11.0  a   =   -9.6 


12.7 


h   =   13.2 


ship  in  rough  seas.   Sea  state 
of  waves  =  40.   Gain  varies.2' 


=  6,  and  relative  direction 


22 


Other  parameters  include  L^=. 51, .50 , . 46 ;  d£=.25,.08; 
e=80.,  a=5.f    $=90.,  6=89.,  h=15. 
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=  12.9(36  =76 
FBRp  =26.9  y 


(b) 

=  30.0 
=  12.1@6=76.0 


(c) 
a    =  60.0 
G    =  9.1(36=76.0 


G    =  12.2(36  =78.0  Gp   =  9.5@6  =79.0 
FBR  =23.4   P       FBR  =  14.2  P 


a    =  90.0 

G    =  3.5(36=76.0 

Gp  =  7.7@ep=i.o 


(e) 


a 
G 
FBR 

6 


90.0 
13.3 
0.0 
180.0 
76.0 


(f) 
a    =  90.0 
G    =  13.3 
FBR  =0.0 
<J>    =  0.0 
6    =  50„ 


Figure  16.   Ten  MHz  3-element  horizontal  array  at 
various  tilt  angles  of  array  axis  (bore  sight). 23 


^Various  antenna  parameters  are:   L^= . 506 , . 500 , . 450; 
d±=.233f  .133;  e=10.,  a=10"3,  4>  =  90w  6=76.  for  (a)  -  (d)  , 


ht=l. 
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Figure  17.   Ten  MHz  3-element 
horizontal  array  operating  over 
land  at  fMH=30«   (For  antenna 
specifications  see  the  ship-ocean 
model  figure  14.) 


fMH=30- 


fMH=35 


Figure  18.   Thirty  MHz  3-element 
horizontal  array  operating  over  land, 
above  and  below  frequency  for  which 
array  was  designed.  (For  specifications, 
see  the  ship-ocean  model  figure  14.) 
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Figure  19.   Ten  MHz  5-element 
horizontal  array  over  land. 
Comparison  of  linear  (inter) 
plot  with  log  plot. 


Gp=14.6  db,  FBR  =6.5  db. 
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Figure  20.   Ten  MHz  3-element 

horizontal  array  over  land. 

Comparison  of  linear  (inter) 

plot  with  log  plot. 

GD=12.2  db,  FBR  =17.7  db. 
v  p 
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24various  parameters  include:  L-=.534 . . 500 , .466 , .466, 
.466;  di=.25, .167, .167 ,.167;  a  =10.,  a  =10"2;  $  =90. ;  6  =77.  ; 
h=l. 

^Parameters  for  this  array  are:   L^=. 550 , . 500, .433 ; 
dj=. 233 , . 133;  all  other  parameters  are  the  same  as  in 
footnote  24. 
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E .   SUMMARY 

There  is  a  somewhat  higher  gain  and  an  increasingly 
higher  front-to-back  ratio  as  the  element  spacing  widens 
in  the  case  of  the  optimal  two  element  horizontal  array 
possessing  a  reflecting  parasite  as  compared  to  the  optimal 
array  with  a  director  parasite. 

The  vertical  multi-element  array  will  produce  a  high 
front-to-back  ratio — with  values  in  the  neighborhood  of 
those  of  the  horizontal  array;  however,  if  resistance  is 
to  be  maintained  above  20fi  it  appears  that  the  resulting 
optimal  vertical  array  will  not  provide  as  high  a  front-to- 
back  ratio  as  the  optimal  horizontal  array.   The  optimal 
design  of  a  multi-element  horizontal  array  provides  both 
a  maximum  front-to-back  ratio  and  a  near  resonant  input 
impedance,  whereas  the  best  design  of  a  multi-element 
vertical  array  represents  a  compromise  between  these  factors, 
e.g.  a  trade-off  sacrifice  exists. 

The  optimal  design  arrived  at  for  a  horizontal  array 
becomes  a  suboptimal  design  when  operated  vertically,  and 
vice  versa. 

Upon  adding  or  eliminating  directors  from  an  array  which 
has  been  designed  for  optimal  performance—with  the  stipula- 
tion that  parameters  of  all  directors  be  kept  uniform,  the 
performance  of  the  array  will  no  longer  be  optimal. 

Operation  of  the  horizontal  array  ©ver  sea  water  environ- 
ment provides  a  higher  gain  with  respect  to  operation  over 
land,  but  over  fresh  water  provides  a  lower  gain  than  over 
land. 
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Gain  is  considerably  higher  in  the  case  of  the  vertical 
array  and  slightly  higher  with  the  horizontal  array  when 
they  are  radiating  from  a  height  of  2\    as  compared  with  IX. 

Gain  varies  considerably  when  the  platform  upon  which 
an  antenna  operates  is  in  heavy  seas. 

The  gain  of  a  horizontal  array  aimed  vertically  into 
the  sky  is  lower  than  when  it  operates  with  no  tilt  angle 
whatsoever. 
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IV.   RECOMMENDATIONS 

At  the  beginning  of  this  study  the  existing  program 
which  had  been  written  for  eight  different  single  element 
antennas  permitted  the  compilation  of,  at  maximum,  only 
three  different  antennas — the  computer  memory  was  exceeded 
when  more  options  were  added.   Now  the  revised  program 
allows  all  nine  antennas  to  be  compiled,  with  the  eight 
unused  antenna  subroutines  residing  externally  on  drum. 
This  permits  consideration  of  further  expanding  the  main 
program  capabilities  inasmuch  as  the  effective  computer 
memory  has  more  capacity. 

In  regard  to  expanding  the  program's  capability  for 
solving  array  problems,  it  is  recommended  that  the  graphic 
input  portion  of  the  program  be  extended  to  allow  for  more 
element  length  and  space  parameters.   This  would  permit 
study  of  long  arrays  (arrays  having  a  cumulative  spacing 
that  total  more  than  one  wavelength) . 

It  is  recommended  that  an  optimization  routine  be  added 
to  the  main  program  which  would  encompass  one  or  more  of 
the  following  criteria:   gain,  front-to-back  ratio, 
resonance  of  the  driven  element  (Z22)  /  thickness  of  the 
element,  and  impedance  (or  standing  wave  ratio) .   For  this 
purpose,  a  more  powerful  computer — which  gives  faster 
solution  time — might  be  preferable. 
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APPENDIX  A 
DESCRIPTION  OF  PROGRAM  OPERATING  PROCEDURES 

The  program  for  the  system  of  antennas  as  originally 
written  has  been  changed  to  accommodate  the  Yagi-Uda. 
Accordingly,  the  program  was  rewritten  in  part  to  obtain 
more  effective  use  of  limited  information  storage  in  the 
computer.   Redundancies  between  antenna  programs  have  been 
reduced,  the  shipboard  simulator  and  the  gain  pattern 
routine  was  rewritten  to  accommodate  the  Yagi-Uda  and  to 
save  storage,  and  the  graphic  text  segment  was  changed  to 
satisfy  the  need  for  additional  parameter  inputs  and  outputs. 
The  program  remains  functional  for  all  antennas  listed  below. 
Figure  21  shows  how  the  parameters  and  the  output  appear  to 
the  operator  at  the  CRT  graphics  display  console. 

To  operate  the  program,  first  ready  the  Graphics  device 
to  be  used  by  following  simple  lab  instructions  for  setting 
up  "Gated. "   The  program  may  be  quickly  loaded  into  the 
computer  by  mounting  the  program  "Dump"  tape  on  tape  unit  #1 
and  loading  a  few  cards — BOOT, A JOB,  AAGT ,  AASSIGN  Xl=MTlA, 
and  ARERUN.   Next,  press  the  usual  operating  buttons  on  the 
9300  console — IDLE,  RESET,  CLEAR  FLAGS,  RUN,  and  CARDS.   The 
cards  and  tape  are  then  read  by  the  computer  and  within  five 
seconds  the  computer  console  will  print  a  request  for  the 
graphics  device  number  that  the  user  desires  to  interface 
with  the  computer.   The  user  then  responds  by  typing  at  the 
IDEV=1*  and  a  carriage  return.   If  a  mistake  is  made  in  typing 
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start  over  again.   If  the  message  was  successfully  typed  but 
for  some  reason  the  text  did  not  appear  on  the  graphics  device, 
the  quick  load  procedure  may  be  used  by  pressing  the  following 
computer  buttons:   IDLE,  RESET,  STEP,  RUN;  then  retype  the 
message  IDEV=1*.   The  graphics  device  should  now  contain  the 
program  text,  providing  all  equipment  is  functional. 

Each  non-zero  parameter  entry  must  consist  of  a  field- 
length  of  fewer  than  four  numbers,  plus  a  decimal,  e.g. 
parameters  should  not  exceed  four  characters  of  information. 
Upon  completing  a  parameter  entry,  press  the  graphics  TTY 
carriage-return  and  the  blinking-light  pointer  will  sequence 
to  the  location  of  the  next  parameter  input.   If  a  parameter 

is  zero  it  need  not  be  typed;  press  the  carriage-return  and 

25 
go  to  the  next  entry  location.     Correction  to  a  typing 

mistake  may  be  made  only  on  the  four-character  field  string 

of  the  parameter  where  the  light-pointer  is  located.   If  the 

error  is  identified  after  sequencing  to  the  next  location  the 

program  may  either  be  executed  as  is  if  no  serious  program 

error  will  result,  of  if  the  error  can  cause  abort  within  the 

SDS-9300,  either  use  the  fast-restart2"  procedure  or  reload 


"In  addition  to  the  need  for  a  decimal  to  be  entered 
with  all  data,  a  blanked-out  space  (one  that  has  been  skipped) 
is  the  same  as  a  zero. 

2°This  procedure  is  similar  to  the  initial  start-up: 
type  the  "Reset  Gated"  messages  at  the  graphics  teletype,  and 
at  the  SDS-9300  press  IDLE,  RESET,  STEP  RUN  and  in  response 
to  the  console  message  listed  by  the  SDS-9300  type  IDEV=1* 
(device  #1) . 
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27 
the  core  "dump"  tape  again.     Either  procedure  takes  less 

than  30  seconds  whereas  problem  solution-time  can  take  longer 


Figure  21.   Photograph  of  graphic 
display  showing  parameter  input/output. 


27  . 

This  procedure  is  identical  to  the  initial  set-up,  in 

that  the  program  dump  tape  is  read  into  the  computer.   Press 

the  button  number  "32"  at  the  computer  console,  and  type 

ARERUN  at  the  console.   This  causes  the  dump  tape  to  rewind 

and  be  read  in  again.   Once  again  the  computer  should  request 

Idev,  and  the  user  should  respond  by  typing  IDEV=1*. 
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Observe  the  parameters  shown  in  the  photograph.  The 
following  is  a  description  of  these  parameters  to  assist 
the  user  in  specifying  his  antenna. 

ANTN:  Antenna 

1.  Tilted  Dipole 

2.  Vertical  Whip 

3.  Vertical  Whip  with  Ground  Screen 

4.  Inverted  L 

5.  Sloping  Long-wire 

6.  Sloping  Vee 

7.  Horizontal  Rhombic 

8.  Vertical  Half-rhombic 

9.  Yagi-Uda 

ELE  LENG:   Element  length  (meters) 

If  a  single  element  antenna  (ANTN  1  through  ANTN  8) 
has  been  specified,  the  first  available  parameter  space 
should  be  used  for  the  element  length  and  the  other  four 
spaces  to  the  right  left  blank.   If  a  Yagi-Uda  (ANTN  9)  has 
been  chosen,  the  following  order  must  be  observed: 
ELE  LENG 


1 

Reflector 

Driven 

Directors 


The  user  may  choose  to  design  the  array  without  a  reflector, 
in  which  case  the  first  location  at  the  left  must  be  left 
blank.   Any  number  of  directors  up  to  three  may  be  selected, 
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however  any  blanked  directors  must  be  contiguous  from  the 
right.   That  is,  skipping  a  director  space  then  entering 
data  into  the  next  space,  will  result  in  the  problem  being 
misinterpreted:   the  first  zero  director  marks  the  end  of 
the  array  and  determines  the  number  of  elements  present. 

HGT/WDTH :   Height  and  width  (meters) 

A  single  element  antenna  requires  only  that  the  height 
be  specified;  the  last  four  spaces  may  be  left  blank,  since 
they  will  be  ignored.   If  the  Yagi-Uda  is  chosen,  select  the 
parameters  according  to  the  following  format: 

HGT/WDTH 


Height    \  \  y^^^ 
Spacings 

The  first  spacing  is  for  the  distance  between  the  reflector 

and  driven  elements.   Subsequent  spacings  are  for  the  distance 

between  respective  directors. 


e 

PHIP:   $  '  Antenna  parameter.   Normally   <t>'=0   will  be  used 
except  under  dynamic  ship  simulation,  where  the  angle  will 
specify  the  array  axis  relative  to  the  ship's  head.   See 
figure  5. 

THEP:   0'  Antenna  parameter.   6'=0   for  a  vertical  YAGI-UDA 
and   9'=90   for  horizontal. 

FREQ:   Frequency  (MHz) 
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EPSL:   Dielectric  constant  of  earth 
SGMA;   Conductivity  of  earth  (mho/m) 
Typical  values  are: 


Sea   i 

Fresh 

Wet 

Dry 

Water  ■ 

Water 

Earth 

Earth 

80   ! 

80 

5-30 

2-5 

3-5    ; 

10' 

"2,  lO"3 

10"1  -  10~3 

10"4,  10~5 

SMUL:   Multiplier 

1.  Multiplies  o    by  . 1 

2.  Multiplies  a  by  .001 

For  numbers   <  3  this  parameter  may  be  used  as  a  multiplier 
of  the  diameter  of  the  array  elements.   The  element  diameter 
is  used  in  determining  the  element  spacing  for  computation 
of  self  impedance.   When  SMUL  >  3  the, expression  used  for 
spacing  is 


VT 


or 


2L  ,  si 


nee   r 


_  *i 


200  200. 

(see  the  computer  program  instruction  that  precedes  #1925) 

PHI:   $   Observation  parameter. 

THET :  e   Observation  parameter. 

Generally  the  observer  wishes  to  position  himself  at 
the  point  of  maximum  gain,  and  take  the  azimuth  and  elevation 
sections  there.   The  first  solution  of  the  problem  will  show 
where  this  point  is  located  by  listing  the  respective  angles 
under  PK:   "Peak".   The  user  may  then  enter  these  values  for 
PHI  and  THET,  and  solve  the  problem  again. 
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LOG;   Log  pattern 

If  left  blank  linear  patterns  will  result.   If  "1." 
is  entered,  a  previously  entered  pattern  via  light  pen 
will  be  erased.   If  "2."  is  entered,  a  log  pattern  will 
result. 

ISTH;   Store  the  horizontal  pattern 

ISTV:   Store  the  vertical  pattern 

Leave  blank  except  to  store  the  pattern  about  to  be 
calculated  for  purposes  of  future  comparison  with  another 
pattern  by  means  of  "double-exposure".   Enter  "1."  to 
store. 

IRCL:   Recall  a  stored  pattern 

Left  blank  until  patterns  have  been  stored  by  a 
previous  problem  solution  using  ISTH/ISTV.   When  comparison 
is  desired  with  the  solution  of  the  problem  being  entered, 
enter  "1".   Ensure  that  if  recall  is  desired  ISTH  and  ISTV 
are  left  blank,  otherwise  comparison  will  be  made  of  the 
solution  with  itself. 

ALPH;   Array-axis  angle 

The  tilt  angle  (or  the  angle  of  elevation  with  the 
ground)  of  the  Yagi-Uda  array  axis.   This  is  the  slope 
angle  with  the  ground  for  the  two  sloping  antennas  (ANTN  5 
and  6)  as  well  as  for  the  vertical  half-rhombic  (ANTN  8) , 
and  the  minor  angle  between  elements  of  the  horizontal 
rhombic  (ANTN  7) . 
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ISEA:   Sea  state 

For  the  tilted  dipole,  vertical  whip,  sloping  long- 
wire,  and  Yagi-Uda.   Left  blank  unless  iterative  solutions 
under  dynamic  sea  conditions  are  desired.   The  program  is 
written  for  a  sinusoidally  varying  sea  with  solutions 
produced  at  thirty  degree  intervals.   The  peak  amplitude 
of  roll  is  8  degrees  per  sea  state,  and  the  peak  pitch  is 
3  degrees  per  sea  state. 

ICRS :   Course  of  the  seas 

Enter  the  direction  of  the  seas  relative  to  the  ship's 
heading . 

GAIN :   Gain  at  the  observation  points  specified  by  <)> ,  6  . 

The  peak  gain,  identified  under  "PK"  is  also  given. 
The  <£  and  e  which  locate  the  peak  gain  are  listed  adjacent 
to  the  values  specified  for  the  problem  (discussed  above 
under  PHI  and  THET) . 

FTBR:   Front-to-back  ratio 

The  ratio  of  the  peak  forward-azimuth  to  peak  reverse- 
azimuth  gains.   (This  is  not  a  comparison  between  the  90 
and  270  degree  gain  values).   See  footnote  3. 
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COMPUTER   PROGRAM 
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